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One of the experimental tests of Lorentz invariance violation is to measure the helicity dependence 
of the propagation velocity of photons originating in distant cosmological obejcts. Using a recent 
determination of the distance of the Gamma-Ray Burst GRB 041219A, for which a high degree of 
polarization is observed in the prompt emission, we are able to improve by 4 orders of magnitude 
the existing constraint on Lorentz invariance violation, arising from the phenomenon of vacuum 
birefringence. 
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I. INTRODUCTION 

On general grounds one expects that the two funda- 
mental theories of contemporary physics, the theory of 
General Relativity and the quantum theory in the form 
of the Standard Model of particle physics, can be unified 
at the Planck energy scale. This unification requires to 
quantize gravity, which leads to very fundamental diffi- 
culties. One is related with the energy of the fundamental 
vacumm state. Another one is the status of Lorentz in- 
variance: the fuzzy nature of space time in quantum grav- 
ity may lead to violations of this fundamental symmetry. 
For the last two decades theoretical studies and exper- 
imental searches of Lorentz Invariance Violation (LIV) 
have received a lot of attention [see e.g. the reviews by 
HH3] . Possible consequences of LIV are energy and helic- 
ity dependent photon propagation velocities. The energy 
dependence can be constrained by recording the arrival 
times of photons of different energies emitted by distant 
objects at approximately the same time [4], e.g. during 
a Gamma-Ray Burst [5] or a flare of an Active Galactic 
Nucleus [BJ. On the other hand, the helicity dependence 
can be constrained by measuring how the polarization di- 
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rcction of an X-ray beam of cosmological origin changes 
as function of energy [7] . 

Recently, an upper limit on the helicity dependence of 
photon propagation has been set using INTEGRAL/SPI 
observations of the polarization of the Crab nebulae [8]. 
In this paper, we derive much stronger constraints from 
a polarization measurement of a Gamma-Ray Bursts 
(GRB), a source at cosmological distance. 

To date only a few polarization measurements are 
available for GRBs. Coburn & Boggs 9| reported a high 
degree of polarization, II = 80% ±20%, for GRB 021206. 
However, successive reanalysis of the same data set could 
not confirm this claim, reporting a degree of polarization 
compatible with zero jTOl [11] . Willis et al. [12] reported 
a strong polarization signal in GRB 930131 (II > 35%) 
and GRB 960924 (II > 50%), but this result could not be 
statistically constrained. GRB 041219A was detected by 
the INTEGRAL Burst Alert System (IB AS; Mereghetti 
IT3"|) . and is the longest and brightest GRB localized by 
INTEGRAL [H] to date Q3]. McGlynn et al. [IB], using 
the data of the INTEGRAL spectrometer [SPI; [T7], re- 
ported a high degree of polarization of the prompt emis- 
sion (II = 68 ± 29%) for the brightest part of this GRB. 
Also, Gotz et al. [TH] have performed a similar measure- 
ment, using this time the INTEGRAL/IBIS telescope, 
reporting variable polarization properties of the burst all 
along its duration. In this paper, we reuse these data in 
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order to measure the polarization in two energy bands 
and check for a shift in the polarization angle as a possi- 
ble effect of Lorentz Invariance Violation. 

In section|TTJ we will present the INTEGRAL/IBIS ob- 
servations of polarization during the burst, and the re- 
cent measure of the burst distance we obtained thanks to 
far infrared observations of its host galaxy. In section |HI| 
we will derive constraints on LIV from these observations 
and we will conclude in section IIVI 

II. OBSERVATIONS OF GRB041219A 

A. INTEGRAL/IBIS observations 

Thanks to its two position sensitive detector layers IS- 
GRI [H] (made of CdTe crystals and sensitive in the 
15-1000 keV energy band), and PICsIT [20^ (made of 
Csl bars and sensitive in the 200 keV-10 MeV energy 
band), IBIS can be used as a Compton polarimeter [21]. 
The concept behind a Compton polarimeter is the po- 
larization dependency of the differential cross section for 
Compton scattering 
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where Tq is the classical electron radius, Eq the energy of 
the incident photon, E' the energy of the scattered pho- 
ton, 9 the scattering angle, and the azimuthal angle 
relative to the polarization direction. Linearly polarized 
photons scatter preferentially perpendicularly to the inci- 
dent polarization vector. Hence by examining the scatter 
angle distribution of the detected photons 

iV(</») = S[l + a o cos(2(0-0o))], (2) 

one can derive the polarization angle PA = 4>q — tt/2 + 
nir and the polarization fraction n = ao/aioo, where aioo 
is the amplitude expected for a 100% polarized source, 
derived by Montecarlo simulations of the instrument. 

To measure the polarization, we followed the same pro- 
cedure described in Forot et al. [22] that allowed to suc- 
cessfully detect a polarized signal from the Crab nebula. 
One important difference, anyway, is that the Crab mea- 
surement was integrated over several observation periods 
and a long time (> lMs), while this measurement here 
integrates over only a few seconds. So we have checked 
if instrumental azimuthal variations, potentially washed 
out by the long Crab observation, could be important in 
the GRB case. To do so, we have simulated ten GRB-like 
observations using data from the INTEGRAL Payload 
Ground Calibrations campaign. This campaign, made 
with the full flight model of the telescope, was done with 
standard radioactive non polarized sources. We use in 
this work the calibrations made with a 113 Sn source with 
a peak energy of 392 keV, the closest to the energy bands 



we are considering here. To simulate GRBs, we divide the 
data set in ten subsets containing each a similar number 
of Compton events as in the observations we show in fig- 
ure[T] In all these subsets, we found systematics at a level 
up to 15% maximum, with a polarization angle around 
40°. We made another tests considering spurious events, 
as described in Forot et al. [22 , and again we found a 
systematic azimuthal variation of less than 9%, with a 
polarization angle of 180°. Both value are far from the 
ones we derived below for the source and give confidence 
that our results are not due to instrumental azimuthal 
variations, even if such variations effectively exist. 

To derive the Gamma-Ray Burst flux as a function of 
<j>, the Compton photons were divided in 6 bins of 30° 
as a function of the azimuthal scattering angle. To im- 
prove the signal-to-noise ratio in each bin, we took ad- 
vantage of the 7r-symmetry of the differential cross sec- 
tion (see Eq. [2]), i.e. the first bin contains the photons 
with 0° < <j) < 30° and 180° < < 210°, etc. The de- 
rived detector images were then deconvolved to obtain 
sky images, where the flux of the source in each bin is 
measured by fitting the instrumental PSF to the source 
peak. We finally fitted the polarigrams to Eq. [2] using a 
least squares technique (see Fig. [2| in order to derive ao 
and (j>o , and the errors on the parameters are dominated 
by the statistics of the data points. 

We computed the scatter angle distribution into two 
energy band in order to detect a possible polarimetric 
angle shift with energy, reminiscent of a possible LIV 
effect. The two bands were chosen to be [200-250 keV] 
and [250-325 keV] where the source has merely the same 
signal to noise ratio. In figure [T] we show the measured 
evolution of the polarimetric angle shift between these 
two energy bands, along the burst duration. These shifts 
are all consistent with zero with a mean value of 21° ±47°. 

Also, as an exemple, we show in figure [2] and figure [3j 
the portion of the GRB light curve where the polarimet- 
ric signal was strong, that is starting at 01:47:02 U.T. 
until 01:47:12 U.T. (P9 interval in Gotz et al. [18]). A 
modulated signal is seen in the two energy bands, corre- 
sponding to n = 55% for the first band, and n = 82% for 
the second ones (see fig{3]for error contours). To evaluate 
the goodness of our fits, we computed the chance proba- 
bility (see Eq. 2 in Forot et al. |2"2"]) that our polarigrams 
are due to an unpolarized signal, and reported these val- 
ues in Fig. [2j The corresponding polarimetric angles 
were PA = 80^8 deg. and PA = 45±|| deg., that is 
consistent at the 2a level (see figure [3]). Propagation of 
errors gives an upper limit of 68° at the 90% level, for a 
possible phase shift. 



B. Distance determination 

Gotz et al. |23j performed deep infra-red imaging of the 
GRB region using the WIRCam instrument at the 3.6 m 
Canadian French Hawaiian Telescope (CFHT) at Mauna 
Kea. Thanks to multi-band (Y JHK S ) imaging they were 
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FIG. 1. Evolution during the burst duration of the polari- 
metric angle shift measured between the [200-250 keV] and in 
the [250-325 keV] energy range. The mean value, 21° ± 47°, 
is consistent with zero. 








50 100 150 

Polarization angle 

FIG. 3. Contour plot of the composed error on polarisa- 
tion angle and polarisation fraction for the two energy bands, 
during the P9 time interval. The X cross shows the best fit 
position for the first energy band. Contour levels are at the 
67, 90, and 95 % level (from white toward black fill). The + 
cross and dotted, dashed, and dash-dotted lines show the best 
fit parameters obtained for the second energy band, consistent 
at 95% with the first energy band best fit. 



200 - 200 keV 




250 - 325 keV 




Angle (degrees) 



III. CONSTRAINTS ON LORENTZ 
INVARIANCE VIOLATION 

On general grounds, Lorentz violating operators of di- 
mension N = n + 2 modify the standard dispersion rela- 
tions E 2 — p 2 + m 2 by terms of the order of f n p n /Mp~[ 2 
where Mpi is the reduced Planck scale (sa 2.4 10 18 GeV), 
used as a reference scale since LIV is expected to arise 
in the quantum regime of gravity. In order to account 
for the severe limits on LIV, one therefore usually only 
considers operators of dimension greater or equal to 5 
which provide corrections which are tamed by at least 
one inverse Plank scale. 



A. dimension 5 operators 



FIG. 2. Scatter angle distribution of GRB04I219A in the 
200-250 keV and in the 250-325 keV energy band during 
the P9 time interval (see text). These distributions give the 
source count rate by azimuthal angle of the Compton scat- 
tering, and are consistent with a highly polarized signal. The 
chance probability of a non-polarized signal is reported in 
each panel. The polarization angles derived from these distri- 
butions are consistent within 68° (see text). 



able to identify the host galaxy of GRB 041219A, and to 
give an lower limit to its photomeric redshift of z = 0.02, 
at the 90% confidence level. This implies a luminosity 
distance of 85 Mpc, assuming standard cosmological pa- 
rameters (fl m = 0.3, £l\ — 0.7, and H = 70 km/s/Mpc). 



If we restrict our attention to pure electrodynamics, 
there is a single term of dimension 5 which gives cor- 
rections of order p 3 /Mpi and is compatible with gauge 
invariance and rotational symmetry |24j : 



L = 
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where n M is a 4-vector that characterizes the preferred 
frame and F 1 *" = \& vpa F prj . The uniqueness of this 
term makes the analysis somewhat model-independent 
(see however below). 

The light dispersion relation is given by (E — hu> and 
p = hk): 
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where the sign of the cubic term is determined by the 
chirality (or circular polarization) of the photons, which 
leads to a rotation of the polarization during the propa- 
gation of linearly polarized photons. This effect is known 
as vacuum birefringence. 

Since we have the approximative relation: 



w± = |pU/l± 
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(5) 



the direction of polarization rotates during propagation 
along a distance d by an angle: 



A0( P ) 



ijJ+(k) — cj_ (fc) 



k 2 d 



2M 



pi 



(6) 



For GRB041219A, if we set A0(k) = 47°, derived from 
the measures we made along the burst duration, and the 
lower limit luminosity distance reported above of d = 85 
Mpc = 2.6 10 26 cm, corresponding to z=0.02, we get an 
upper limit on the vacuum birefringence effect: 



2M Pl A6(k) 
{k\ - k\)d 



1.1 10" 



-14 



(7) 



B. dimension 6 operators 

Although the dimension 5 operator (|3| is unique, it is 
physicaly relevant only if there does not appear operators 
of lower dimensions. In the general case however, radia- 
tive corrections induced by dimension 5 operators lead 
to dimension 3 operators through quadratic divergences 
which contribute an extra factor A 2 , where A is an ultra- 
violet cut-off of order Mpi . Barring extreme fine tuning, 
one needs a symmetry argument to cancel such terms. 
Supersymmetry appears to be the only symmetry which 
cancels such contributions [35] . It is true that supersym- 
metry is broken in nature but, in the case where it is 
softly broken, quadratic divergences contribute a factor 
M|, where Ms, the scale of supersymmetry breaking, is 
of the order of a TeV 2 , thus providing an extra factor 
(Ms /Mpi) 2 ~ 10~ 30 compared to untamed dimension 3 
operators. 

Unfortunately, the dimension 5 operator of Q is not 
compatible with supersymmetry [25) . We therefore have 
to resort in this case to dimension 6 operators. We thus 
assume that the light dispersion relation is given by: 
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We can derive then the approximate relation: 
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which implies that the direction of polarization rotates 
during propagation of: 



A0(k) = 



w+(fc) — w_ (k) 



k 3 d 



(10) 



Again, for GRB041219A, if we set A0(k) = 47° and 
the luminosity distance d — 85 Mpc = 2.6 1 26 cm, corre- 
sponding to z=0.02, we get an upper limit on the vacuum 
birefringence effect: 



2M Pl A0(k) 
(kl - k\)d 



2.6 10 8 



(11) 



This is still too large an upper bound to be really con- 
straining since one expects couplings at most of order 
one. 



IV. CONCLUSIONS 

Using a recent determination of the distance of 
GRB041219A [53] for which a high degree of polarization 
is observed [TH1 US] ; we were able to increase by 4 orders 
of magnitude the existing constraint on Lorentz invari- 
ance violations, arising from birefringence. Turned into 
a constraint on the coupling £ of dimension 5 Lorentz vi- 
olating interactions, that is of corrections of order k/Mpi 
to the dispersion relations, this gives the very stringent 
constraint £ < 10~ 14 . Most presumably, this means that 
such operators are vanishing, which might point towards 
a symmetry such as supersymmetry in action. In that 
case, the pressure is on the next corrections of order 
(k/Mpi) 2 corresponding to operators of dimension 6. We 
showed that, although astrophysical constraints are not 
yet really constraining, they are getting closer to the rel- 
evant regime (£ of order 1 or smaller) . Our result can be 
compared to the limits derived using the possible energy 
dependence of the group velocity of photons in distant 
GRBs derived with Fermi [5]. However, written in the 
same way as we did in this work, this limit was £ < 0.8 
that is much less constraining. 
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